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BUCK!, TNG COF A LONG SQUARE TUBE IN TORSION 
, AND COMPRESSION 

By Bernard Budl snaky, Manuel Stein, and Arthur C. Gilbert 


SUMMARY 


The buckling of an infinitely long square tube under combined 
torsion and compression is investigated by means of an enact energy 
method utilizing Lagrangian multipliers. An interaction curve is 
obtained from which it is possible to determine the amount of one 
loading required to produce buckling when a given amount of the other 
loading is present. 


INTRODUCTION 


The local buckling of a long, thin-^wall, square tube subjected to 
a combination of torsion and longitudinal ■ compression is investigated 
theoretically in the present paper. The walls of the tube are considered 
to be infinitely long flat plates continuous over ncudeflecting line 
supports at the comers. An exact theoretical analysis of the problem 
by means of the lagrangian multiplier method, presented in detail in 
the appendix, is used to derive an interaction curve which giveB the 
combinations of torsion and compression required to produce local 
buckling. 
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width of tube wall 


wall thickness 
Young’s modulus 
Poisson’s ratio 


plate stiffness in bending 
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compressive stress 
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shear stress 


torque 


/ b^t'N 

compressive— stress coefficient \o J 
shear— stress coefficient t — — I 


ratio of shear stress present to critical shear stress for 
pure torsion 

ratio of compressive stress present to critical compressive 
stress for pure compression 

plate coordinate parallel to length 

plate coordinate parallel to width 

deflection normal, to plane of plate 


a^bjj^Cjjijdni Fourier coefficients 
a*,!)* Lagrangian multipliers 

X half wave length 


internal handing energy' 

external work of compressive stress 


external work of shear stress 


integers 


EESUITS AND DISCUSSION 


The critical combination of c ompr essive and shear stress of a 
long square tube (Bee fig. l) is given by the formulas 


a = kc 


r = k s 
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where r is related to the total torque T by the formula 



An interaction curve for the critical combinations of the stress 
coefficients k c and k 0 is presented in figure 2; the theoretical 
analysis from -which this curve is obtained is given in the appendix. 

Essentially the interaction curve is a curvilinear portion Joined 
to a straight portion; it is of interest to compare this curve with the 
interaction curve for an isolated simply supported plate alBo shown in 
figure 2 (see reference l). The disparity in behavior of the two curves 
may be explained in terms of nodal patterns. 3h the case of the isolated 
plate, the inclination as well as the spacing of the transverse nodal 
lines adjust themselves so as to cause the buckling stress to be a 
minimum. Thus, as relatively more compression is applied, the nodal lines 
become less inclined, until, for the case of pure compression, they are 
straight and perpendicular to the edges of the plate. However, in the 
case of the square tube, the inclination and spacing of the nodal lines 
are always constrained to be such that a nodal line b9 continuous all 
the way around the tube. But, as relatively more compression is applied, 
at a certain ratio of compressive to shear stress, the nodal lines 
suddenly cease to be inclined and became straight and perpendicular to 
the tube comers (and thus still remain continuous around the tube). 
Buckling with this type of nodal pattern corresponds to the straight 
portion of the interaction curve. The closeness of the two. curves in 
the shear-predominating range is due to the fact that the nodal pattern 
for the isolated simply supported plate in pure shear happens to be one 
that would very nearly be continuous around a square tube made up of 
four such plates. ^ 

The interaction curve for the square tube is shown in stress— ratio 
form in figure 3 and is compared with a parabolic interaction curve that 
is shown in reference 1 to hold very closely for an isolated plate 
having equal elastic restraint of any magnitude along the edges. The 
comparison shows clearly that, at least in the compression— predominating 
range, it would be unduly conservative to consider th.9 walls of the tube 
to behave as isolated elastically restrained plates. 

CONCEDDHTG BEMAEKS 


A theoretically computed interaction curve for the buckling of an 
infinitely long square tube in torsion and compression is presented. 
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From, this curve it can "be concluded that an appreciable amount of torsion 
may "be present without in any way reducing the compression required for 
buckling. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aer on autics 
Langley Field, Ya., September 10, 1$A8 
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THEORETICAL ANALYSIS 


For the purpose of the present analysis, the square tube is 
idealized into four infinitely long flat plates continuous over non- 
deflecting line supports. Evidently two distinct types of buckle 
pattern must be considered. The first type (see sketch) given by 


w = sin ~ sin ^ 
b b 


represents the case in which the nodal lines are straight and do not 
advance longitu dinally as they proceed around the tube: 



-k- b -*| 


For this buckling mode, the tube will buckle only if 


a = k 


b^t 


This type of buckling is represented by the straight portion of the 
interaction curve (see figs. 2 and 3)j lor this type of buckling, the 
shear stresses do no external work. 
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In the second type of ‘buckling deformation, the nodal lines are 
inclined and advance longitudinally exactly one full wave length as they 
traverse the four vails of the tube an shown in the following sketch: 


7 



X -*| 


From the general remarks concerning Fourier series in appendix B of 
reference 2, it is seen that in the region (0,X), (0,2b) the deflec- 
tion may he given by 
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sln^l 
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sin 
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2b 




ut^l, 3 ^5^ 1 


b cos 
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+ cos 
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m=l,3j5j • 


nwy 

8in 2b” + 



jft=l,3j5.> • • • 


dm 008 



(1) 


(Single, rather than double, Fourier series are required because the 
deflection iB Binusoidal along any line in the infinite direction. ) 

This function satisfies the requirement that 

v(x,y) = -v(x,y+2b) 

Additio nal conditions that must be satisfied are 

v(x,y) = v(X-x,b-y) (2) 

and 

v(x,y) = -^w(X-x,3b-y) 


(3) 
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Equations (2) and (3) are, respectively* conditions uf symmetry in 
the region (O^X), (O^b) and antisymmetry in the region ( 0 ,X)j (h J 2 h). 

mi 

2 

These conditions are fulfilled "by making b_ = a^— l) and 

ffidL 

2 

= — c (-i l) . Thus the deflection function becomes 


w = sin S a m 


nt=lj3j5^ • • ■ 


sin + (-1) 


mr-1 

2 


cos 


mity 
2b . 


+ cos 2S. ^ 


y ®= 1 , 3 , 5 ,' 


Cjn. 


mr-1 

cos BE 


( 4 ) 


This function will be used in an exact stability analysis by the 
Lagrangian multiplier method as described in reference 2 . 

The bo undar y conditions of zero deflection along the idealized 
supports (comers) 


w(XjO) = w(x,b) = 0 


lead to the constraining relationships 


00 mr -1 

XI %(-!> 2 - 0 ( 5 ) 

nt=lj 3^5j • • • 


and 


00 m — 1 

X- <^(-1) = 0 

nt=lj 3^ 5 } • • • 


(6) 
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The Internal energy Y and the external work of the stresses T c 
and T 0 are given by the expressions 



Substituting the Fourier expansion of v (equation (4)) into the 
energy expressions yields 


jb y- 

32b 2 p 3 **=1,3,5,... 


(4 + m 2 p 2 ) 2 (a m 2 + Cjjj 2 ) 


(7) 


and 


^ 51 (=» 2 + ° m 2 ) 

2b p m=l,3,5,... 






nt=l,3,5, • 


maj^CmC-l) 




(8) 


(9) 
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■where 


P = F 


k c = a 


t 2 t 

k s = T -gr 
it D 


The function to he minimized is 


F = T — T c — T 0 — a* 


m?=l,,3,5,. 


(“!) % - n 1 




Minimizing equation (10 ) -with, respect to the coefficients aj and Cj 
yields 


p ^ 2 

= (4 + J 2 P 2 ) a, - l6k 0 2 a.» + let-P 3 ^-!) c, - a(-l) = 0 




(k + J 2 p 2 ) 2 cj - l6k c 0 2 Cj + l6k B P 3 j(-l) 2 aj - T) ( 1 ) 2 = 0 (12) 
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where 


3 - • • • 


a = a* 


l6t 2 p 3 


and 


i) = V 


l6t 2 p 3 


After simplification., equations (ll) and (12) 10001119 

3-1 

A 3 a 3 + B 3 C 3 = a( - l} 

and 

3-1 

B 3 a 3 + a 3 c 3 = 1 ^ -1 ) 

where 


Aj = (4 + J 2 p 2 ) 2 - l6k c p 2 


and 


M, 

Bj = l6k B p 3 3(-l) 
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Then. 



xfcl 

2 2 
Ajtt(-l) -B^riC-l) 



( 13 ) 


and 




2 2 
A t jT](--l) - B ^aC-l) 





Substituting far a^ and bj In the constraining relationships ( 5 ) 
and (6) yields 



nt=lj3j5 j • • • 


Ajjft, — Bmjj 

-Am 2 ~ Bm 2 


( 15 ) 


and 



21=1, 3^5 j • • • 


■Am 1 ! ®m a ' 
% 2 "^2 


( 16 ) 


Tor the lagrangian multipliers, a and t|, to have values other than zero, 
the condition that must be satisfied is 


2l 


Am 


iut=l,3,5, ... Affl^ 





( 17 ) 



12 


NACA TO Ho. 1751 


Equation ( 17 ) can "be factored, into 


J22 : 



m— 1 


1 

+ ®m 


= 0 


( 18 ) 


and 


£ 


m=l,3,5, ... Am ®m 


0 


(19) 


Both equations (l8) and (19 ) will give identical results. The posit iTe 
or negative sign merely indicates the direction of shear. 

Equation (l8) niay he written 


OO 

0 . n 

“= 1 ^ 3 , 5 , • • • 


1 

ffi -1 

(4 + Ei2p2) 2 _ l63£cp 2 + Igkgmp^-l) 2 


( 20 ) 


The critical combination of compressive— stress coefficient k c and 
shear-stress coefficient k s for an infinitely long square tube can he 
calculated from equation (20 ). For a specified value of one of the 
stress coefficients and wave-length ratio P, the corresponding value of 
the other stress coefficient that will satisfy equation (20) can he 
obtained. This procedure is used for several different values of P 
until a minimum value of the corresponding stress coefficient is obtained. 
The curvilinear portion of interaction curve was drawn using these 
minimum values. (See table 1.) 
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